Electron spin dynamics and electron spin resonance in graphene 
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A theory of spin relaxation in graphene including intrinsic, Bychkov-Rashba, and ripple spin-orbit coupling 
is presented. We find from spin relaxation data by Tombros et al. [Nature 448, 571 (2007).] that intrinsic 
spin-orbit coupling dominates over other contributions with a coupling constant of 3.7 meV. Although it is 1-3 
orders of magnitude larger than those obtained from first principles, we show that comparable values are found 
for other honeycomb systems, MgB2 and LiCg; the latter is studied herein by electron spin resonance (ESR). We 
predict that spin coherence is longer preserved for spins perpendicular to the graphene plane, which is beneficial 
for spintronics. We identify experimental conditions when bulk ESR is realizable on graphene. 
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Introduction. The discovery of graphene JJl] stimulated 
enormous interest due its fundamentally and technologically 
important properties. One potential application is in spintron- 
ics iQl, i e. when the electron spin degree of freedom is uti- 
lized as information carrier The principal parameter govern- 
ing spintronic usability is the spin relaxation time (also re- 
ferred to as spin-lattice relaxation time), Ts, which character- 
izes how an injected non-thermal equilibrium spin state de- 
cays. For realistic applications, Tj longer than 10-100 ns is re- 
quired. A general, often cited concept is that "pure materials 
made of light elements" can reach this limit. The huge mo- 
bility of charge carriers in graphene (approaching 10^ cm^A' s 
(jstl), the light nature of carbon, and the low-dimensionality 
of this material are the reasons for the high expectations for 
its spintronic applications. This is supported by the long spin 
relaxation time in light metals such as e.g. Li |01 or in low- 
dimensional conductors [5]. 

Therefore it came as a surprise that Ts as short as 60-150 ps 
are observed in spin transport experiments on graphene iaLZl], 
which renders it unusable for such applications. The under- 
standing of this experimental result is therefore of great im- 
portance. Theory of spin relaxation are split into two different 
classes: materials with inversion symmetry (e.g. Na or Si) 
and to materials where the inversion symmetry is broken ei- 
ther in the bulk (e.g. III-V semiconductors such as GaAs) or 
in two-dimensional heterostructures. The Elliott- Yafet (FY) 
theory ^ explains the former case, where only intrinsic 
(i.e. atomic) spin-orbit coupUng (SOC) is present, L^, and 
predicts that spin (Fj = ft/Tj and momentum relaxation rates 
(F — H/t, t is the momentum relaxation time) are propor- 

tional: F. = cti-^T. Here ai = 1..10 is band structure de- 
pendent ij], A is the energy separation of a neighboring and 
the conduction band. 

The relaxation for broken inversion symmetry is explained 
by the Dyakonov-Perel (DyP) theory. It applies either when 
the symmetry breaking is in the bulk, (the Dresselhaus SOC 
ifioll . Ljx) or when it happens for a heterolayer structure (the 
Bychkov-Rashba SOC lullill, Lbr)- The DyP theory shows 
that the spin and momentum relaxation rates are inversely pro- 



portional: Fs = aD/BRio/BR/r. where od/br « 1- 

A hnk between the FY and the DyP was found recently 
II13II : for metals with inversion symmetry but rapid momen- 
tum scattering, the generalization of the EY theory leads to 

Fs = cti A^i+r^r, which gives a DyP like spin relaxation when 
F > A. 

Three sources of SOC are present in graphene: intrinsic, 
BR type (due to the symmetry breaking by a perpendicular 
electric field), and the ripple related (which is due to the in- 
evitable ripples in graphene). However, the role and magni- 
tude of these SOC parameters is a debated issue. Estimates 
for the intrinsic SOC ranges two orders of magnitude; 0.9- 
200 ^eV [14-16], whereas value of the BR SOC appears to 
be settled to 10-36 ^eV per V/nm (Refs. jllj] and llj, re- 
spectively). The effect of the substrate for the spin relaxation 
is also unsettled Given this debate, a description is re- 
quired which enables comparison with the spin transport data. 

Here, we present the theory of spin relaxation in graphene 
including intrinsic, BR, and ripple spin-orbit coupling. We an- 
alyze the spin transport data from Refs. |l6l[7lll8t] and we find 
that the intrinsic SOC dominates the relaxation with a large, 
unexpected magnitude. We discuss two similar honeycomb 
systems; MgB2 and LiCe, and show that they exhibit similar 
intrinsic SOC. The result predicts a strong anisotropy of the 
spin relaxation time. We study the feasibility of bulk electron 
spin resonance (ESR) spectroscopy on graphene and pinpoint 
experimental conditions when it is possible. ESR would allow 
a direct, spectroscopic measurement of Ts (Ref. |fl9|]), which 
underlines its importance i20|] . 

Experimental. We prepared Li intercalated HOPG graphite 
by the "immersion into molten Li" method 112 ill . The golden 
color of the samples attested the LiCe intercalation level li22ll . 
Freshly cleaved samples were sealed under He in quartz tubes 
for the ESR experiment. 

Spin relaxation in graphene. Low energy excitations 
around the K point of the Brillouin zone are described by a 
two-dimensional Dirac equation: 



H = VF{axPx + (yyPy), 



(1) 
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with the v-p w 10^ m/s Fermi velocity iU]. The spin-orbit 
interaction in graphene is given by fl^: 



LbR + -^ripple (r) 



{axSy — (JySx), (2) 



where Li, Xbr, and Lrippie are the SOC's of the intrinsic, BR, 
and ripple terms, respectively, irippie(r) is Gaussian corre- 
lated random variable, (Lripple(r) -^ripple (r')) ^ — I"')- 

The spin relaxation rates induced by these 
SOC's are additive in lowest order provided 

max(ii,LBR,irippie) <max(r,^): 



Ts - Ts.i + Ts. 



BR 



Ts, ripple 



(3) 



The contributions from the intrinsic (F^ i), BR (F^ br), and 
ripple (Fs.rippie) relaxation rates are obtained using the Mori- 
Kawasaki formula similar to that used in Ref. ll23ll considering 
the conical band structure and the K, K' degeneracy: 



Lf arctan(/^/F)^ 
27r/z • F) 



{26^, !_ + £br 
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I' =11, or ± is the spin polarization direction with respect to 
the graphene plane; e.g. i/ =|| in the spin transport experi- 
ments [6]. Here, /i is the chemical potential and F) = 

In (^^^-^T-^ + ImI ^ f arctan (^^]) is the pseudo 
chemical potential (I? « 3 eV is the cutoff in the continuum 
theory) which appears in the expression of the density of states 
(DOS), p{ii, F), with finite ^ and F: 



2AeM(/i, r) 



(7) 



with Ac = 5.24 A^/(2 atoms) being the elementary cell and 
F) is measured in units of states/eV • atom. 
The intrinsic contribution disappears when spins are polar- 
ized perpendicular to the plane and the BR and ripple terms 
have a 2:1 anisotropy for the ±:|| directions. For the intrinsic 

part, Fs.i « (2/7)2 ^ when /i ^ F, which is an Elliott- Yafet like 
result with a, — 1 since the band-band separation, A = 2/i. 
In the vicinity of the Dirac point, DP, (i.e. /i « and F finite) 

it returns a Dyakonov-Perel like result of Fs.i = 4i„( jyp) r- 
This is in agreement with the generalized Elliott- Yafet theory 
which predicts a similar crossover when the momentum scat- 
tering rate overcomes other energy scales [13]. Interestingly, 
the intrinsic contribution can be well fitted with a Lorentzian: 
Ts.i « a'j^^, where a' « 0.2..0.4 and F'/F w 1..2 for 
typical values of ^ and F. 
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FIG. 1: a) Experimental (symbols, from Ref. fW]) and calculated 
spin lattice relaxation rates, Fs, as a function of /i in graphene for in- 
plane spin polarization. Upper(lower) solid and dotted curves are the 
maximal(minimal) estimates for the intrinsic and BR contributions 
with SOC values from Refs. fl4l - fl^ . respectively. Dashed curve 
is the ripple contribution with Lrippic from Ref. I14n . The upmost 
thin solid line is a fit to the data as explained in the text, b) The 
same experimental data shown as Ts along with the fit (solid curve). 
For comparison, Ts calculated with F = 12 meV (dashed curve) is 
shown. Arrows depict the crossover of the DyP and EY mechanisms 
as a function of p. 



The BR term is only present if a perpendicular electric field, 
E, is applied, which induces a BR SOC of Lbr = kE with k 
values between 10 ifTsIl and 36 /ieV/(V/nm) (l^. The electric 
field changes /i through: /i = ^JjmWv^ where n = f3E is 
the carrier density and /3 = 0.22 (V • nm)^^ for Si02 gate 
insulator [24]. This yields the BR SOC as a function of /x: 



iBR(A*) « • 3.4^^7^. 

The ripple relaxation contribution depends on F only if 
/i ^ F, where it resembles an EY relaxation: Fj ^ppie cx 
tippler ln(i?/F). 

Analysis of the spin transport data. In the following, we an- 
alyze the available spin transport data |0, 0, [Tsl] in the frame- 
work of the above calculation. Values of Ts = 60..125ps 
were found around the charge neutrality point (depending 
on the sample), with a typical F w 75meV llST . Fig. [T] 
shows the measured and calculated spin relaxation rate data 
for =||. F = 75meV, that is independent of /i, was 
used for the calculated curves. First principles calculations 
of the intrinsic SOC scatter more than two orders of magni- 
tude with values of 0.9 /xeV H S S), 24 /ieV fil, and 
200 ^eV Values for the BR SOC, Xbr = kE, vary be- 
tween K = 10.. 36 /ieV/(V/nm). This gives rise to the minimal 
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and maximal estimates for both types of the contributions as 
shown in Fig. [T] The ripple SOC was estimated to be 17 /ieV 
inRef. Q. 

Clearly, the first principles based relaxation rates fall short 
of explaining the experimental observation. Of the three con- 
tributions, only the intrinsic one has a /i dependence that mim- 
ics the experiment, whereas the other two shows the opposite. 
It may appear that a fit to the data is ill defined, given the rela- 
tively large number of free parameters (F and 3 L's). However 
to our surprise, the fit consistently yields the same, robust set 
of parameters, irrespective of starting values or the method 
used (least squares fitting or combined with a simulated an- 
nealing), which are: Li = 3.7(1) meV, Lbr = irippie = 0, 
F = 120(5) meV. This robustness originates from the qual- 
itative difference between the /i dependence of the different 
contributions. The obtained values satisfy the criterion for the 
perturbative approach and the value of F determined herein is 
in agreement with that obtained in Ref. [18]. 

The intrinsic SOC opens a bandgap of Li in the excitation 
spectrum [ 15,. 16.1 therefore it is natural to ask: why is not this 
gap observed experimentally? Two interrelated answers are 
in order: first, best quality samples to date are ballistic only 
on the (sub)micron scale, giving a momentum scattering rate 
of the order of meV's (or bigger), which can mask the gap 
11271]. Second, charge inhomogeneities (the so-called puddles) 
prevent us from reaching the Dirac point, the average minimal 
charge density is estimated 02 811 as 10^ cm~^, which gives an 
average /i ~ 4 meV, capable of overwhelming the obtained 
gap- 

The present analysis allows for the design of graphene 
based spintronic devices. For spins polarized perpendicular to 
the graphene plane, the intrinsic contribution vanishes thus re- 
sulting in a substantially longer spin relaxation time. For spins 
polarized in the graphene plane. Fig. [T] shows that around 
the Dirac point purer samples (i.e. smaller F) decreases Tj 
rather than increasing it, thus deteriorating performance. This, 
somewhat counterintuitive phenomenon, is the consequence 
of the Dyakonov-Perel like behavior of the intrinsic contribu- 
tion around the DR 

The large value obtained for the intrinsic SOC is surprising 
as it is an order of magnitude larger than the largest theoretical 
estimat e Ml 611 a nd up to 3 orders of magnitude larger than other 
results il4l25ll26ll ). However, given that the experimental fi 
dependence of Fs dictates the dominant role of the intrinsic 
coupling, Li yields necessarily a large value. In the following, 
we consider two similar systems, MgB2 and Li doped graphite 
and show that therein similar values of the intrinsic coupling 
are obtained. 

In MgB2, the boron atoms form a honeycomb lattice with 
four p-shell electrons, such as in graphene, which highlights 
the similarity of the two materials. Therein, an intrinsic SOC 
of Li(MgB2) = 2.8 meV of the tt orbitals was found i23lLI t 
was shown by Griineis and coworkers [2^ and confirmed 130(1 
that alkali atom intercalated graphite is an excellent model 
system of biased graphene as the two-dimensional electron 
dispersion is retained due to the weak interlayer coupling. The 
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FIG. 2: High temperature ESR linewidth in HOPG LiCs (full sym- 
bols) and linear fit to the data (solid line). We show similar data from 
Ref. (31I1 (open symbols) on a LiCe powder sample. Inset shows the 
schematics of the LiCe band structure according to the PES measure- 
ments 1.29,, 30.1 and the A = 1.65 eV parameter. 




Li intercalated stage I graphite compound LiCe 12211 is par- 
ticularly suitable to determine the intrinsic SOC as Li is the 
lightest alkali metal and its contribution to the spin relaxation 
is undetectable 101 • 

In Fig. ID, we show the temperature dependent ESR 
linewidth, AB, for an HOPG LiCe along with previous data 
on a powder LiCg sample |31] and schematics of the band 
sti-ucture. A linear fit to the data yields AB = 0.205 [mT] -f 
T X 6 • 10^^ [mT/K]. Of these terms, the temperature de- 
pendent one is associated with the homogeneous broaden- 
ing, ABhom, due to SOC, which gives Fj — .g/iBA_Bhom 
(g ~ 2 is the g-factor, /iB is the Bohr magneton) and is 
Fs = 2.1 ■ 10"^ eV at 300 K. Since Li doped graphite re- 
sembles biased graphene i29[ [30ll . the above theory of the 

intrinsic SOC applies, i.e. Fs = ^VT. With the values of 
A = 1.65 eV 130] and a typical F(300K) = 4.4meV [22], we 
obtain ii(LiC6) = 1.1 me V. Although it is debated whether 
SOC in graphite is applicable for graphene the similar 
result for these three systems leads us to conclude that the in- 
trinsic SOC is properly determined in graphene. 

Detectability of ESR on graphene. With the SOC param- 
eters and the theory of spin relaxation at hand, we assess 
the feasibility of ESR spectroscopy on graphene. It is de- 
termined by the sample amount, the magnitude of the spin- 
susceptibility, and the ESR linewidth. The ESR signal is pro- 
portional to the amount of magnetic moments: xoVB/^q 
where B is the magnetic field, xo is the volume spin- 
susceptibility (dimensionless in SI units), V is the sample vol- 
ume, and iiQ is the permeability of vacuum. For graphene with 
area A, the amount of magnetic moments is xo.gr^B / fiQ with 
Xo,gr having a unit of meters. The Pauli spin-susceptibility 
of graphene is xo.gi- = MoMbPIMi T)^, N is the number of 
carbons and the DOS, p{fi, F), is given above. 

ESR spectrometer performance is given by the limit-of- 
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FIG. 3: Limit of ESR detection for graphene as a function of /i and 
r in units of the graphene area (upper panel) for an in-plane mag- 
netic field. Arrows show the maximum chemical potential by gate 
bias and by chemical doping and solid curve indicate the 100 mm^ 
area border Expected ESR linewidth, Ai?|| (lower panel), solid lines 
show two selected linewidths, 1 and 10 mT. 



detection (LODq) i.e. the number of S = 1/2 non-interacting 
spins at 300 K which give a signal-to-noise of = 10 

ratio for AB = 0.1 mT linewidth, and 1 sec/spectrum-point 
time constant. For state-of-the-art spectrometers LODq = 
10^° spins/0.1 mT. The spin-susceptibility of such spins is 
Xcurie — fJ'OfJ'B^ / i^^T), wherc spins occupy a volume 
of V, which gives an LOD for graphene: 



LOD„ = LODq 



/(AS) 



26meV x piii,T) 



(8) 



in units of number of carbons. Here, 26 meV is the thermal 
energy at 300 K and the /(AB = 0.1 mT) = 1 function de- 
scribes that the ESR S/N decreases as 1 /AB if AB < 1 mT 
(the typical magnetic field modulation limit) and as 1/ [ABY 
if AB > 1 mT. For n > T, the DOS is well approximated by 
p(/i) = 0.0385/z • [states/eV^atom] which yields a compact 
result: LODgr w 1000 LODq • f{AB)/n with ^ in eV units. 

Clearly, a sizeable DOS and narrow linewidth are prereq- 
uisites to observe ESR on graphene. Large DOS can be 
achieved by moving /i away from the DP or by inducing de- 
fects. The latter yields, however, increased scattering thus 
larger linewidth. Shifting /i by a gate bias is limited to 
~ 0.2 eV due to breakdown in the most common Si02 insu- 
lator around iJ w 0.1 V/nm. With chemical doping using K, 
•apto^^ 1.35 eV can be achieved ll29ll . The ESR linewidth 



is expected to be strongly anisotropic with a minimum, AS^, 
for a perpendicular magnetic field, its magnitude however re- 
mains unknown. In Fig. |3], we show the calculated LOD for 
graphene as a function of /i and F in units of the graphene 
area along with the calculated linewidth for an in-plane mag- 
netic field. The LOD can be two orders of magnitude smaller 
for a perpendicular magnetic field if the corresponding ESR 
linewidth is an order of magnitude smaller. Therefore ESR 
experiments should be attempted with the perpendicular ori- 
entation first. This experiment would yield directly the magni- 
tude of the BR and ripple relaxation contributions from AB±. 

An important benchmark, that indeed the intrinsic ESR sig- 
nal of graphene is observed, is an angular dependence of the 
ESR linewidth: AB{9) = sin^ (61) AB|| + cos^{0)AB^ as 
a function of the asimuth angle, 9. Finally, we note that 
the anisotropy could reconcile the narrow ESR linewidth in 
the perpendicular geometry [i20il with the short in the spin 
transport experiment 0] . 

In conclusion, we presented a theory of spin relaxation 
in graphene which takes into account intrinsic, Bychkov- 
Rashba, and ripple spin-orbit coupling induced spin relax- 
ation. Analysis of spin relaxation data show that the intrinsic 
contribution dominates the relaxation with a coupling constant 
that is orders of magnitude larger than theoretical estimates 
but it is not unusually large compared to other honeycomb 
systems. The result predicts a large anisotropy of the spin re- 
laxation. We presented under what circumstances bulk ESR 
spectroscopy can be observed in graphene. 
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